Unraveling the fluctuations of animal motor activity 
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Human motor activities are known to exhibit scale-free long-term correlated fluctuations over a 
wide range of timescales, from few to thousands of seconds. The fundamental processes originating 
such fractal-like behavior are not yet understood. To untangle the most significant features of these 
fluctuations, in this work the problem is oversimplified by studying a much simpler system: the 
spontaneous motion of rodents, recorded during several days. The analysis of the animal motion 
reveals a robust scaling comparable with the results previously reported in humans. It is shown 
that the most relevant features of the experimental results can be replicated by the statistics of the 
activation-threshold model proposed in another context by Davidsen and Schuster. 

PACS numbers: 87.19.L, 89.75.-k , 89.75.Da 
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I. INTRODUCTION 

The timing between consecutive human actions, or 
even spacing the most simple and inconsequential mo- 
tor activities [J, [3, is known to be scale-free, and de- 
spite recent efforts their mechanisms remain poorly 
understood. The interest is further underscored by its 
potential usefulness as statistical markers for the detec- 
tion and follow up of human neurological disorders. The 
lack of plausible models to account for these statistical 
features calls for alternatives able to identify the essen- 
tial underlying mechanisms. To that end, it may be ad- 
vantageous to uncouple the interference of memory and 
other cognitive factors, present in the experiments with 
human subjects, by analyzing a more elementary pro- 
cess, namely the records of long-term spontaneous motor 
activity of laboratory animals. 

In the present work, the activity of rats was recorded 
during several days, and the experimental data series 
analyzed from the perspective of a point process. The 
results show striking similarities with the scaling demon- 
strated to be exhibited by humans in more elaborated 
settings. The paper is organized as follows: In the next 
section the experimental details are described. Section 
III contains the statistical analysis, first for the inter- 
event times, then for the rates of motion events and fi- 
nally for the variance of counts. Section IV describes an 
activation-threshold model which is able to replicate the 
most relevant experimental observations. Finally, section 
V closes the paper with a discussion of further implica- 
tions of the present results. 



II. EXPERIMENTAL SETTING AND DATA 
RECORDING 

Male 4-month-old Wistar rats were kept in a sound- 
proof room temperature (~ 20°C— 22°C) and humid- 
ity (~ 50 — 80 percent) conditioned chambers. Ani- 
mals were individually isolated into transparent cages of 



25cmx 25cm X 12cm with food pellets and tap water ad 
libitum. They were exposed to a 24-hour cycle of light- 
dark conditions: cycles of 12 hs of light (fluorescent lamps 
with intensity of 300 lux) and 12 hs of darkness (dim red 
light with intensity less than 0.1 lux) [5]. 

Animal's movement was monitored with infrared 
activity-meters, able to scan each rodent housing cage 
and to report, at a frequency of 1 Hz, animal movements 
even in absence of locomotion. Given the sensor's sen- 
sitivity the detection includes even minute head motion, 
grooming, etc. An example of the recording for the first 
three days is plotted in Fig. [1] where changes from activ- 
ity to immobility, and viceversa, are marked by a vertical 
bar, with a resolution of 1 s. In the bottom plot of Fig. [2 
we also show the group average of the activity rate R, 
which is a coarse graining of the raw data, exhibiting the 
well known circadian rhythmicity. 
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FIG. 1: Three days activity plot obtained from six labo- 
ratory animals ("A" trough "F") exposed to a cycle of 12 
hours of light (6AM-6PM) and darkness (6PM-6AM), indi- 
vidually housed and continuously monitored by an infrared 
device scanning the field at a rate of 1 Hz. The bottom time- 
series, depicts the group average activity (R), computed with 
a binwidth of 1 min and normalized. 
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Our interest here focuses exclusively on understanding 
the dynamics of the irregular fluctuations and not on the 
circadian periodicity. 



III. STATISTICAL PROPERTIES 

The spatial and temporal resolution of our experimen- 
tal setting is such that, from the original (binary) data 
recorded, we can precisely estimate the duration of mo- 
tion (sequences of I's) and immobility (sequences of O's) 
episodes. Thus, the animal motor activity can be inter- 
preted as a point process where the beginning of a mo- 
tion event occurs at a definite discrete time ti . The point 
process can be then specified by the sequence of event 
times ti or, alternatively, by the series of increments (or 
inter-event intervals) Ti = ti+i — ti. A typical plot of the 
T-series is presented in Fig. [2] It is evident that, beyond 
the circadian rhythmicity, clustering or bursts of activity 
occur. 
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FIG. 2: Inter-event time n as a function of the event number 
i. (Data from animal A.) 

In order to characterize this point process we will first 
consider the r-series, obtaining its distribution and corre- 
lation structure. Next, rates of events will be inspected. 
Finally, the statistics of event's counts Q as a function 
of observation times will be described. 



A. Inter-event intervals 

Fig. [3] shows the estimation of the probability den- 
sity P{t) through the relative frequency of occurrence 
of inter-event times. It was computed over the nine-day 
activity (binary) recordings, for each one of the six ani- 
mals. We have overimposed the plots from all the animals 
to demonstrate the robustness of the results. It can be 
clearly seen that, from few seconds to several thousands 
seconds (about 1 hour), the distribution of inter-event 
times decays as a power-law (with exponent ^ falling 
within the interval 1.75 ± 0.5 for all six animals). 
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FIG. 3: Normalized distributions of: inter-event times r 
(empty symbols), duration of motion episodes (gray symbols) 
and duration of immobility periods (black symbols), com- 
puted from nine days of continuous recording. Symbols joined 
by dotted lines correspond to the results from each one of the 
six animals. The solid lines correspond to a double exponen- 
tial fit (with characteristic times of the order of 1 s and 4 
s). The dashed line, drawn for comparison, has slope -1.75. 
Inset: representation of the distribution of motion episodes in 
log-linear scale. Statistics computed individually for each of 
the six animals and plotted overimposed 



In contrast, the distribution of duration of motion 
episodes do possess a characteristic timescale (see also 
Fig. [3]). It can be described by a superposition of two 
exponentials with characteristic times of the order of 1 s 
and 4 s, close to the smaller data resolution and to the 
average duration of motion episodes, respectively. Let 
us note that, for human (arm) motor activity data, in- 
stead of two exponentials, a stretched exponential fit was 
reported 0. 

Meanwhile, quiescence intervals are also power-law dis- 
tributed with exponent fi. Because motion intervals are 
in average much shorter than quiet ones, then, the statis- 
tics of inter-event times is mainly dominated by that 
of immobility periods, in particular, sharing the same 
power-law decay. Hence, discrepancies between both his- 
tograms are evident for small time intervals only (see 
Fig. [3]). For comparison, let us note that, for the power- 
law exponent of the cumulative histogram of inactivity 
periods of human (arm) motor activity, the values 0.92 
(control) and 0.74 (depressed individuals), with about 
10% relative error, have been reported 3- 

Identical distributions to those exhibited in Fig. [3] 
were obtained when day and night fiuctuations were an- 
alyzed separately, only differing in the cutoffs, occurring 
at longer activity (inactivity) intervals, at day (night). 
Therefore the fundamental statistical features are inde- 
pendent of the level of activity. 

To assess correlations, first we computed the spectral 
density of r, displayed in Fig. [51 showing that inter-event 
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intervals are not independent. Therefore, the stochastic 
point process of motor activity can not be considered a 
renewal one. Because of the relatively small value of the 
exponent of the r power spectrum, we computed also 
that of the increments, which helps to confirm that it is 
not white noise. Although not shown, we verified also 
that motion and quiescent intervals are anticorrelated. 
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FIG. 4: Long-range correlations of inter-event times. Log-log 
plot of the power spectra for the six series of inter-event inter- 
vals (open symbols) and of their increments (filled symbols). 
Data were logarithmically binned. Here / is in units of the 
inverse event number. Dashed lines are a guide to the eyes, 
with slopes —0.3 and 1.7. Notice that, a = 0.3 and /3 — 1.7, 
verify a = 2 — /3. (Data from six animals.) 



B. Local rates 

By dividing the whole observation time interval in 
W (nonoverlapping) uniform windows of length T and 
counting the number of events, 7V„, in each time window 
n, one obtains the series of counts. The cumulative num- 
ber of events, Af{t) = Yln=i vs. time t = Xl^^i 
is shown in Fig. [5l^a). One observes that the local rate 
(local slope) is not constant as in standard Poisson pro- 
cesses. It varies on time, following not only the slowly 
circadian trend but also other rapidly changing ones. 

In order to depict more precisely rate inhomogeneities, 
we estimated local rates as being the ratio of Nn/T. 
In the scale of half-day, one observes two characteristic 
mean rates associated to the half-periods of low and high 
activity levels. However, shorter time windows T reveal 
a more complex structure of the distribution of rates, 
that is not simply bimodal but has the shape displayed 
in Fig[5jb). It follows as a power-law (with exponent 
7 ~ 0.75) with exponential cut-off. The higher proba- 
bility of small rates leads to longer intervals of immobil- 
ity. Notice also that there is a range of window lengths 
for which the distribution of rates remains invariant, a 
scaling feature typical of a fractal-like process. Besides 
this scaling property, rates fluctuate in time, as exhibited 
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FIG. 5: (a) Cumulative number of events as a function of 
time. A zoom is displayed in the inset, (b) Distribution of 
local rates R computed over non-overlapping time windows of 
lengths T indicated on the figure. The solid line corresponds 
to a power-law with exponent 7 = 0.75, drawn for compari- 
son. Upper inset: the same plot in linear scale. Lower inset: 
timeseries of local rates R for two of the time window lengths 
considered. (Data from animal A.) 



in the lower inset of Fig. [3{b). Nonhomogeneous Pois- 
son processes with (uncorrelated) stochastic rates have 
been considered to explain the emergence of scaling in 
the statistics of inter-event times (see for instance [7]). 
Since in the present case the inter-event intervals are not 
independent, such inhomogeneous Poisson processes can 
be excluded as responsible for this dynamics. 

To estimate rate linear correlations, we performed a 
spectral analysis of the time series of increments /„ = 
Rn — Rn-i- Fig. [Hlis a log-log plot of its power spectrum 
S{f). It is approximately linear on a wide range of biolog- 
ically relevant temporal scales implying that S{f) ~ f^^ 
with /3 ~ 1. Since /3 > 0, consecutive values of the pro- 
cess / are negatively correlated, meaning that increases 
in activity are, on the average, more likely to be followed 
by decreases and viceversa. Shuffling the timeseries of 
increments yields white noise. Recalling that R is the 
integration of /, then the spectral density of the original 
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R timeseries decays as 1//" with a = 2 — f3 
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FIG. 6: Long-range anticorrelations in tiie increments /„ of 
activity rates _R„ computed over 1-min windows. Log-log plot 
of the power spectra for the six / timeseries (symbols) . The 
dashed line with slope /3 = 1 indicates the correlations ex- 
pected for pink noise. (Data from six animals). 



C. Variance of counts 

Common measures for detecting correlations in se- 
quences of counts are the Fano {FF) and Allan {AF) 
factors [gI]. The former is the ratio of the variance to 
the mean of the number of events in each time-window, 
FF = (([7Vfe]2) - (7Vfc)2)/(iVfc), an index of dispersion of 
counts. The latter quantifies the discrepancy of counts 
between consecutive windows, being AF = {{[Nk+i — 
Nk]'^))/{2{Nk)). Although related, they reflect differ- 
ent features, therefore we kept track of both of them. 
In Fig. [71 the two factors are plotted as a function of 
the length of the counting time window T. There is 
a range where they increase as T'^. The exponent of 
FF is d & [0.65,0.73] (that is d ~ fJ- - I, within er- 
ror bars), for the six laboratory animals, while that of 
AF is about 0.1 higher. The discrepancy may be due 
to the fact that FF{AF) tends to the power law from 
above(below). Their power-law behaviors indicate that 
the temporal support of the spikes of activity is a fractal- 
like set with dimension d. Shuffling the series of inter- 
event intervals modifies both factors, mostly by reducing 
the upper bound of the power-law scaling region. This 
indicates that the scaling properties are partially due to 
the distribution of inter-event intervals itself, but also 
that some correlations are associated with the specific 
ordering of the intervals. 



IV. MODELING THE DYNAMICS 

The observed super-Poissonian behavior discussed in 
the previous Section points in the direction of a mul- 
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FIG. 7: Fano and Allan factors as a function of the counting 
time- window length T. Solid lines are the results of fits giving 
slopes d ~ 0.7 and d ~ 0.8, respectively. For comparison, the 
same analysis over the series of shuffled time intervals is also 
shown. (Data from animal A.) 



tiplicative or clustering process, where the occurrence of 
an event increases the probability of a subsequent one 0] . 
A class of clustering Poisson process was introduced by 
Griineis [IQ] . where there is a primary Poisson process 
triggering the occurrence of a sequences of events (clus- 
ters), each following a secondary Poisson process. In each 
cluster the number m of events is a random variable. 

The statistical properties of this two-stage cascade are 
determined also by the cluster size distribution p{m) . A 
special case of interest in the context of 1//" fluctuations 
is p{m) (X for m < A'o, and null otherwise. For 

3/2<77<3, in the limit A'^o — * oo, it was reported 
that the exponent of the variance/mean curve is d ~ 
(7 — 277) /4 while the exponent of the spectral density is 
3 — ry. If r/ = 2, 1// noise is obtained and d = 3/4, which 
is in good accord with the present outcomes. Although 
Griineis et al. clustering Poisson process can reproduce 
some of the observed scalings, it is not clear how such 
a cascade process would precisely be originated in the 
present context. 

The same difficulty applies to other fractal or fractal- 
rate stochastic point processes [1]. Recall also that, for 
many processes cited in Ref . Q , solely the distribution of 
rates is scale-free, while inter-event times present a more 
trivial statistics. In the present case, however, inter-event 
times are not only scale-free but also present long-range 
correlations, indicating non renewal processes. 

It is plausible that animal activity is triggered when 
an internal dynamical state variable reaches some value. 
Without being too specific, an animal can move to eat 
when sugar level reaches some low value, for instance. 
Of course, biological reality will indicate that nothing in 
the judgement of the state variable nor in the thresh- 
old value can be very precise. Therefore, one can imag- 
ine a quantity relaxing towards a fluctuating threshold 
that resets upon crossing it. Variants of this scenario 
are very common in the the literature where fluctua- 
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FIG. 8: Relaxation toward a fluctuating threshold model, (a) 
Plots as a function of time of: the threshold B (black line), 
that fluctuates diffusively (with coefflcient D) between the 
levels and V\ (dashed lines) ; the dynamical variable V (thin 
red line) that describes the relaxation from the excited level 
V2 (dotted line) down to the threshold and is characterized by 
amplitude K and exponent k. (b)-(d) Results of simulations 
performed up to time lO'' in units that correspond approx. 
to 1 s of real time. Model parameters are (Vi, V2, D, K, k) — 
(27,30,0.04,3.0,0.25). (b) Series of inter-event intervals, (c) 
Distribution of inter-event times. Dotted line indicates ^ = 
1.75. (d) Fano and Allan factors vs. counting window length 
T. Dotted line corresponds to d = 0.7. 



tions are introduced in either the threshold level or in 
the activation function [ill, [Hj [H, [l3l- In particular, 
we examine here the variant introduced by Davidsen and 
Schuster [31 j where the threshold fluctuates following a 
(bounded) Wiener process with diffusion constant D. In 
Fig. [Hl^a) , we present a representative example of the fluc- 
tuating threshold & and the relaxation dynamical vari- 
able V as functions of time. 

Fig-[Hb)-(d) show the model dynamics. The distribu- 
tion of inter-event times depends on the specific shape of 
the decay of V{t) , the relaxation process. If the decay has 
the shape V{t) = V2— K{t—tiastY , where tiast is the time 



of the previous adjacent trigger, then the power-law dis- 
tribution of inter-event times has exponent fi = 2 — k [l^ . 
In particular k = 0.25, yields fi = 1.75 close to the ob- 
served values, and as soon as k approaches zero (abrupt 
decay), the exponent can increase up to a value slightly 
smaller than 2. This relaxation ruled by k might explain 
the observation in Ref. [2| of ^ ~ 1.7 in depressed indi- 
viduals and ~ 1.9 in control ones. 

Besides the inter-event time distributions, the other 
main correlation features seen in the data are repro- 
duced as displayed by the FF and AF measures (notice, 
though, the relatively larger dip for short T in the AF). 
It is known that the model dynamics yields 1// spec- 
tral properties [l3|. We have verified (not shown) that 
the inter-event times are power-law correlated and that 
the distribution of rates is scale free, as observed for real 
data, although the exponents are different for the chosen 
set of parameter values. 

Overall, this simple model is able to reproduce the ob- 
served phenomenology and it is known to be robust un- 
der the addition of noise over the activation signal [l3| ■ 
Finally, the duration of individual events could also be 
incorporated easily into the model by an additional inte- 
gration process consistent with the observed statistics. 



V. FINAL OBSERVATIONS 

In summary, the results show that animal motion is 
scale free in all of the relevant statistics analyzed and that 
the scaling is introduced by the length of the inactivity 
pauses as well as its specific ordering. 

These results are robust for all animals studied, and in- 
variant when day and night data are analyzed separately. 
While motion episodes do possess a definite timescale 
and are basically exponentially distributed, the distri- 
bution of inter-event intervals are scale-free over many 
timescales, from a few up to thousands of seconds. Scal- 
ing is also demonstrated in the the rates of occurrence, 
rejecting a non-homogeneus Poisson process as responsi- 
ble for the dynamics. 

At even longer time scales than the ones considered 
here, many human actions exhibit also bursty patterns, 
in which short time intervals of intensive activity are (fre- 
quently) separated by long gaps without activity [4|. In 
that case, the power-law tails of inter-event distributions 
have been attributed to queuing processes (where the 
queue represents the tasks on a priority list). Depending 
on the dynamics of the priority list, different classes of 
universality can emerge |^] . Probably this is not the case 
in the present experiments, where more likely, an ele- 
mentary physiological and/or environmental mechanism 
is behind. 

Some of the scale-free properties here observed, have 
been reported before for human motor activity data as 
well [ij, [^. The quantitative aspects for the scaling 
laws were reported to remain unchanged under usual 
daily activities or periodic scheduled work. These ob- 
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servations were consistent with earlier results from the 
study of timeseries of heartbeat intervals [1, [l^, [lB| from 
healthy humans. Then, it have prompted the sugges- 
tion of multi-scale physiological mechanisms as responsi- 
ble for the presence of long-term correlated dynamics 1] . 
Other related approaches [13] consider that animal mo- 
tion could be visualized as the output of a large nonlin- 
ear dynamical system (e.g., the brain-body-environment 
ensemble) whose repertoire includes the kind of dynam- 
ics observed in these experiments. Finally, the possibility 
that the complexity of the environment in itself influences 
animal behavior needs to be carefully considered [18i] . 

The specific mechanism discussed here, involving a 
fluctuating-threshold controlled dynamics is biologically 



plausible, and not only can replicate observed features 
of the statistical structure of rat motor activity, but with 
appropriate modifications can be used to provide insights 
about the long-term rhythms alterations observed on in- 
dividuals with mood disorders, depression and other neu- 
rological disorders, such as chronic pain. 
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